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The relation between the vector-isoscalar and vector-isovector interactions in the two-flavor Nambu-Jona- 
Lasinio (NJL) model is investigated under the different constraints from QCD. We demonstrate that the flavor- 
mixing can be induced by the mismatched vector-isoscalar and vector-isovector interactions at finite baryon 
chemical potential u and isospin chemical potential uz. The effect of this non-anomaly flavor-mixing on the 
possible separate chiral transitions at nonzero uy is studied under the assumption of the effective restoration 
of U(1)4 symmetry. We find that for the weak isospin asymmetry, the two separate phase boundaries found 
previously can be converted into one only if the vector-isovector coupling gY is significantly stronger than the 


vector-isoscalar one g$ without the axial anomaly. 
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I. INTRODUCTION 


Plotting the map of QCD phase diagram at finite tempera- 
ture T and quark chemical potential u has attracted growing 
interests. Especially, a chiral critical endpoint is predicted by 
some model studies, which separates the crossover and the 
first order phase transition. Currently, the lattice QCD com- 
putations at finite density still struggle in a limited range of 
H [1]. So the existence and location of the critical endpoint 
are still under debate due to the lack of reliable theoretical 
methods for the non-perturbative dense QCD. Theoretically, 
fluctuations of conserved charges, such as net-baryon, net- 
charge and net-strangeness, are predicted to be sensitive to 
the correlation length of the system [2—4] and directly con- 
nected to some susceptibilities [5, 6]. Thus, the experi- 
mental data related to these quantities can serve as powerful 
tools to probe the critical endpoint in heavy-ion collisions: 
the search for such a point is ongoing at RHIC (BES) [7-9] 
and will be performed in the future facilities in GSI (FAIR) 
and JINR (NICA). Moreover, unconventional multiple chiral 
critical endpoints are also proposed: it is found in [10, 11] 
that the finite isospin chemical potential u; may lead to two 
critical endpoints; when considering the color superconduc- 
tivity (CSC), the low-temperature critical endpoint(s) may 
appear due to the interplay between the chiral and diquark 
condensates [12-17]. 

Besides T and u, the U(1)4 anomaly may impact the QCD 
phase transition significantly [18]. This point has been con- 
firmed in model studies or Ginzberg-Landau analyses, where 
the U(1)4 anomaly is usually incorporated by introducing 
the Kabayashi-Maskawa-’t Hooft (KMT) interaction [19-21]. 
The KMT interaction explicitly breaks the U(1)4 symmetry 
and gives rise to the flavor-mixing among light quarks: the 
dynamical mass of u quark may contain contributions from 
both d and s quark condensates; the diquark condensate for 
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u-d pairing at moderate or high baryon density may make 
contribution to the s quark mass. Consequently, the U(1), 
anomaly may affect not only the properties of the traditional 
critical endpoint [22], but also that of the unconventional one: 
the separation of chiral transition due to finite uy [10, 11] may 
be removed by the anomaly flavor-mixing [23]; a new critical 
endpoint at low temperature could be induced in the presence 
of the color flavor locking (CFL) CSC [13]. 


However, the recent lattice calculations indicate that the 
U(1)4 symmetry may be restored obviously near and above 
T. for zero u (24, 25]. The effective restoration of the U(1), 
symmetry would influence the universality class and criti- 
cal properties of the chiral transition [26]. Phenomenologi- 
cally, the model studies suggest that the properties of the con- 
ventional critical endpoint are quite sensitive to the degree 
of U(1), restoration [22]. Moreover, if the anomaly related 
flavor-mixing is very weak near the phase boundary, the two 
critical endpoints due to the isospin asymmetry [10, 11] could 
still be possible due to the decouple of light quarks. 


On the other hand, it is also probable that the non-anomaly 
flavor-mixing can be induced by other ingredients of QCD, 
especially under some condition. The main purpose of this 
work is to study the possible non-anomaly flavor-mixing of 
light quarks and its effect on the chiral phase transition under 
the isospin asymmetry. In particular, we will focus on the fate 
of the two critical endpoints due to the separate chiral transi- 
tions found in [10, 11] with the assumption of the effective 
restoration of U(1)4 symmetry near the phase boundary. 


Our starting point is the four-quark vector interactions with 
different coupling strengths in the isovector and isoscalar 
channels. The effect of vector interactions on the chiral tran- 
sition has been extensively studied in the NJL-type model of 
QCD. A well-known result is that the chiral transition at finite 
4 is weakened by the vector-isoscalar interaction (Wy phy 
and the critical point disappears for strong g5 [12, 15, 22, 27]. 
It is also found that the g% in a proper range can lead to the 
new low-temperature critical endpoint(s) [12] when consid- 
ering the two-flavor CSC, especially with the constraint of 
electric-charge neutrality [15, 17]. In general, g5 and the 
isovector coupling g} are independent coupling constants un- 
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der the chiral symmetry. In the literature, the non-anomaly 
flavor-mixing due to the mismatched vector interactions in 
the vacuum has been discussed in Ref. [28] based on a three- 
flavor NJL model. We extend such a work to finite T and u 
to study the effect of the non-anomaly flavor-mixing on the 
chiral transition with isospin asymmetry. 


Il. EXTENDED NJL-TYPE MODEL WITH 
MISMATCHED VECTOR INTERACTIONS 


A. The general four-quark interaction model with 
mismatched vector interactions under the chiral symmetry 


We start with the following Lagrangian of four-quark inter- 
action model for two-flavor QCD 


P ae (1) 
with 
3 
Lm = 801 > [OYP + Grtiyswy] 
a=0 
3 


= 82 AU yh? + ie ys] 


a=0 
= gale yh) + Tyu ysy] 
= gal yy — (br yyyswyl, (2) 


and 


LY = gsoldet{WC — ys)W] + h.c.} 
= gly — PPY — diyswy + Griy], O) 


where 7°, and 7 refer to the unit matrix and Pauli matrices 
in the flavor space, respectively. The former term Ly in 
(1) is the general Fierz-invariant form of the four-quark inter- 
actions in color-singlet channels which respecting the global 
flavor symmetries of SU(2)y@S U(2)4®U()yeu(1), [29]. 
The latter one pi is the KMT interaction induced by the 
gauge configurations of instanton and anti-instanton [19, 20], 
which only possesses the $ U(2)y@S U(2)4@U(1)y global fla- 
vor symmetries. 

As mentioned, we will focus on the flavor-mixing arising 
from the mismatched vector interactions at finite density. We 
see that three of four independent coupling constants in Le 
are related to the vector and axial vector interactions. Gen- 
erally, the nozero sum g,3 + gy4 implies that the vector cou- 
pling strength in the isovector channel is different from that 
in the isoscalar one. Similarly, the non-vanishing 8,3 — 8v4 
indicates the axial-vector interactions are also mismatched in 
the isovector and isoscalar channels. How the vector cou- 
pling difference gives rise to the non-anomaly flavor-mixing 
at finite density will be detailed in next section. 

The axial-vector interaction may be responsible for the de- 
viation of the chiral magnetic effect in the recent lattice cal- 
culation compared to the analytic formula, as proposed in 
Ref. [30]. Here we mainly study the chiral phase transition in 
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the MFA, the axial-vector interactions in Lagrangian (1) will 
be ignored. Thus we only consider the the following effective 
Lagrangian, 

3 


LR = gar Y OY + eiysWy] 


a=0 
+ gY — (baw)? — (hiyswy + Yriysy)] 
= Sy why - bey wh”, (4) 


where the independent coupling constants are reduced to four. 


B. Unequal vector coupling constants in the mean field 
Hartree-Fork approximation 


Here we stress that the vector coupling difference in the 
MFA can also arise from a very popular version of the NJL 
model [29], 


3 
LY = ga J [Oy + Oriy] 
a=0 
+ga? - GP? - Wiysh)” + Utiyswy’] 
3 
= gv Y [OEV + Or Yuy], (5) 
a=0 
in which only one vector coupling g, is adopted. In the 
Hartree approximation, there is no difference between the 
coupling strengths of the two vector interactions at the mean 
field level for Lagrangian (5). 

However, the effective vector couplings (in the sense of 
direct interaction) in the isoscalar and isovector channels will 
differ from each other if the Fock contribution is also consid- 
ered. For a four-fermion interaction, the Fock contribution 
can be easily evaluated according to its Fierz transformation 
[29]. Taking into account the exchange terms, the effective 
direct four-quark interactions of the Lagrangian (5) take the 
following form 


4 4 4 
L£ ) sff-direct = L! ) + Dag 


= (Bs + 80+ SU) + GiystU)"] 
+ (831 ~ 8:0 EP + Wish) 


3 


= gv DIG yy? + Or ys] 


a=0 


= ( &v 4. Tags 

Ne 2Ne 
sl 
Ne 
where N, is the color number of the quarks. The effective 
Lagrangian (6) clearly shows that the exchange terms give 
rise to the vector coupling difference in the Hartree-Fock ap- 
proximation (HFA), which is at the order of O(1/N.) com- 
pared to g,,; and g,. Note that the similar result in a three- 
flavor NJL model has been given in [28]. 


uy, 


aliii 
S a Oog, (6) 
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If both g, and gs; in (5) originate from the color current- 
current interaction (yap), they fulfill the relation g, = 
gsı/2 according to the Fierz transformation. In this case, the 
vector coupling difference shown in (6) becomes 


s v 2 §s1 
gv = 8-8 = wie) = Wo (7) 


This equation indicates that g% in (4) may be larger than gy 
and their difference is considerable compared to g, or gs; for 
N: = 3. 


C. Constraints on the vector interactions from the lattice 
chiral curvatures 


Even Eq. (7) implies that the coupling gY is weaker than g5, 
it is also possible that g} may be stronger than g5. This can be 
understood from the curvature difference for the chiral phase 
transition at finite baryon and isospin chemical potentials ob- 
tained in recent lattice calculations [31]. 

For small baryon and isospin densities, the chemical 
potential dependence of the pseudo-critical temperature for 
the chiral crossover can be expressed as 


Tc(Ug, Hi) = Te + Agi + Bip; + OMG Ha), (8) 


where Te is the chiral pseudo-critical temperature at zero 
quark chemical potential (In this subsection, 44 and yu; are 
used to refer to the quark baryon and isospin chemical poten- 
tials, respectively). Notice that T.(ug, 4i) is an even function 
of uqi [32]. So at the order of Hoe we can expand Teu?) as 


2 
Kofi 
Telg) = Tell = kqi yr} (9) 
C 
where the two chiral curvatures are defined as 


ict (10) 


The lattice QCD simulation in [31] suggests that the curva- 
ture x, is about 10% greater than kj. 

Recently, the role of g} on the determination of x, has been 
studied in a Polyakov-loop enhanced three-flavor NJL model 
[22]. It is found that x, decreases with g} and to reproduce the 
lattice x, the g} must keep relatively larger value compared to 
gy. The authors of Ref. [22] then propose the lattice x, can be 
used as a useful constraint on g5. 

We can directly extend this idea to determine x; by replac- 
ing 4q with u;. As will be demonstrated in the next section, 
the coupling g/ influences the curvature «x; in the similar way 
as gp does on ky. In particular, x; and x, obtained at the 
MFA of the two-flavor NJL model will take the same value 
for g) = g5. In other words, the lattice curvature difference 
between x; and x, can be regarded as an useful evidence for 
the unequal vector coupling strengths. 

Since the two-flavor lattice calculation in [31] indicates that 
Ki is less than x,, we thus infer that g} may be larger than 
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g; near the chiral phase boundary for zero and small quark 
chemical potential. Following the spirit of Ref. [22], our nu- 
merical study suggests that g} is about 10% larger than g% near 
T, according to the lattice curvatures in [31]. Note that this 
conclusion is quite different from the estimation given in (7). 


D. Constraints on the vector interactions from the couplings 
of vector mesons to nucleons and lattice susceptibilities 


In Ref. [33], it is argued that the ratio of the couplings of 
w and p mesons to nucleons can be used as a constraint on 
the vector coupling difference. In the chirally broken phase, 
the empirical value for this ratio is given by gunn/8pnn ~ 3, 
whereas in the chirally symmetric phase it is expected to be 
one. It is then proposed that the ratio g/g% is located in the 
range from 1/3 to 1. 

In addition, another quite similar estimation is given in 
Ref. [34], where the vector coupling difference is expressed 
as the function of two susceptibilities y, and yz under some 
assumptions. Using the lattice data for these susceptibilities 
as input, it is found that g) is always less than g5: their dif- 
ference is quite large below Te which approaches zero rapidly 
above T. for zero chemical potential. 

All the arguments given in the above subsections suggest 
that the vector interactions are repulsive (namely, g} and gý 
are all positive), but the relation between the gf and gy re- 
mains uncertain. In the following study, g5 and g} in (4) will 
be treated as the free parameters due to these uncertainties. 


II. VECTOR-INTERACTION INDUCED 
FLAVOR-MIXING AND THE THERMAL DYNAMICAL 
POTENTIAL AT FINITE BARYON AND ISOSPIN 
CHEMICAL POTENTIALS 


In this section, we shall demonstrate that the vector 
coupling difference can lead to non-anomaly flavor-mixing at 
finite baryon and isospin densities. 

The full Lagrangian of two-flavor NJL model with the 
interaction (4) reads 


L = Wid," + yojt— to) y + LO. (11) 


where the quark chemical potentials are introduced and ño = 
diag(m,, mq) is the current quark matrix. We shall adopt the 
isospin symmetric quark masses with m, = mg = mo. The 
fin the Lagrangian (11) is the matrix of the quark chemical 
potentials which takes form 


A Hu HHI 
2 = ; 12 
g ( Ha ) ( H+H ) pal 
with 
Hu + Ha HB Hu — Ha ôu 
= z= — d = = —. 13 
H 7 z @d p 5 7 (13) 


In (13), 4g (ur) is the baryon (isospin) chemical potential. 
Note that the definition of the isospin chemical potential in 
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the quark level is different from that in the nucleon level. For 
more details on the role of isospin symmetry energy in nu- 
clear matter, the reader can refer to [35-38] and references 
therein. 

At finite densities, the quark chemical potentials are shifted 
by the vector interactions. Here we use w’ to denote the modi- 
fied quark chemical potential. Note that the u quark density is 
different from the d quark one under the isospin asymmetry. 
The shifted quark chemical potentials take the form 


Hid) = Hud) ~ 28)(Pu + Pa) — 28;(Oud) — Pa) 
= Hud) 7 L(g; T 8 )Puld) E (g; a 8 Pdi)» (14) 


or 


K =u- 285(Pu F Pa), Hy = ur 28) (Pu = Pa) (15) 


where 


Pud = YY) (16) 


is the u (d) quark number density. Eq. (14) clearly shows 
that due to the vector coupling difference, not only p, but 
also p4 give contribution to the effective chemical potential of 
u quark, and vise versa. This implies that the flavor-mixing 
arises due to the vector interaction. As mentioned, this mix- 
ing has nothing to do with the axial anomaly. The modified 
chemical potentials can also be rearranged as Eq. (15), which 
indicates that u and u; are shifted by the isoscalar and isovec- 
tor vector interactions, respectively. 

Formally, the non-anomaly flavor-mixing shown in (14) 
for the modified chemical potentials is quite similar to 
the anomaly flavor-mixing for the constituent quark masses 
induced by the instantons, namely: 


Mua) = Mo — 4251 Pua) a 48 .0bau)s (17) 
where 
ud) = OCCA (18) 


is the u (d) quark condensate. 
Using the conventional technique, the mean field thermal 
dynamical potential of the Lagrangian (11) is expressed as 


Q(T, Hu, Ma) = 
X, Q(T, u p Mp) + 28O + $2) + 4B 2buba 
f=u,d 
= (85 + 8O + P3) — 2(85 - 8Y)Pupa» (19) 


where Qo(T, u p; My) is the contribution of a quasi-particle 
gas of the flavor f which takes the form 


Q(T, u ps Mp) = 


ap ; 
-onet | SP [mi + exp(-Ey— a! pD] 


+ Inf + exp(-(Ey + A/T] 


ap i 2 2 
-2N; f ays ENA -P (20) 
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with the quasi-particle energy Ey = pt M4. The A in 
Eq. (20) is the parameter of three-momentum cutoff in the 
NJL model. We see that besides the modified chemical po- 
tential Hp» the flavor-mixing due to the vector coupling differ- 
ence is also explicitly demonstrated in Eq. (19) via the direct 
coupling between p, and pg. 

Minimizing the thermal dynamical potential Eq. (19), the 
motion equations for the mean fields ¢,, a, Pu and pg are 
determined through the coupled equations 


Q Q Q Q 
eg ag Mug “sg. (21) 
Opu Opa 


This set of equations is then solved for the fields ¢,, 
Éd, Pu and pa as functions of the temperature and chemical 
potentials. 


IV. FATE OF THE SEPARATE CHIRAL TRANSITIONS 
WITH NON-ANOMALY FLAVOR-MIXING 


As mentioned, the separate chiral transitions due to finite 
Hi (10, 11] can be removed by the flavor-mixing induced by 
the axial anomaly [23]. Since the instanton density may be 
suppressed significantly near the phase boundary, we revisit 
this problem by taking into account the non-anomaly flavor- 
mixing due to the mismatched vector interactions. 

For comparison, we follow the notations in Ref. [23] and 
introduce two parameters œ and g, which are defined as 


851 = (d — Q)85, 852 = A8s, (22) 


here œ means the ratio of the KMT interaction in the scalar- 
pseudoscalar channel, which is treated as a free parameter 
in the following calculations. The other model parameters, 
namely the current quark mass mo, the scalar coupling con- 
stant g, and the three-momentum cutoff A are all adopted 
from [23]. 


A. Fate of separate chiral transitions under the weak isospin 
asymmetry without the axial anomaly 


The role of the mismatched vector interactions on the sepa- 
ration of chiral transition at finite T-u under the weak isospin 
asymmetry is investigated by switching off the KMT interac- 
tion. We focus on whether the two critical endpoints found 
previously could be ruled out by the non-anomaly flavor- 
mixing without the help of the axial anomaly. 

We first study the cases for gf > g5 with a fixed small 
coupling gẹ = 0.2g, under the weak isospin asymmetry 
ôu = —20 MeV (The typical value of du in heavy ion colli- 
sions may be within this range, as estimated in [23]) 3. The 
T-z phase diagrams for varied g} are shown in Fig. 1. For 


3 Note that the jy defined in [23] corresponds to the ôu in our notations 
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80 : : 
(a) 
60+ g= 0.29, 
gy = 0.49; 
D 40L H =-20MeV 
= 
H 
20 
(0) A fi 1 . 
300 350 400 450 
u (MeV) 
80 T T 
(b) 
60 g = 0.2g; 4 
g% = 0.6g; 
© 40L Sp =-20MeV J 
= 
m 
0 
300 350 400 450 
H (MeV) 
80 
t (c) 
60+ ON = 0.29, 4 
| gy = 0.68g, 


4o- du =-20MeV 4 


300 350 400 450 
u (MeV) 


Fig. 1. (Color online) The T-u phase diagrams for varied vector- 
isovector coupling g} at du = —20 MeV without the axial anomaly. 
The vector-isoscalar coupling g is fixed as 0.2g,. The solid line 
stands for the first-order chiral boundary. 


gi = 2.0g%, Fig. 1(a) shows two separate first-order phase 
boundaries, which correspond to the chiral transitions for the 
u and d quarks, respectively. For g; = 0.6g,5, Fig. 1(b) shows 
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that only one first-order chiral boundary emerges at the low 
temperature, but it splits into two lines at the relatively higher 
temperature. So there are still two critical endpoints. Fur- 
ther increasing g) to 0.68g,, Fig. 1(c) displays that only one 
phase boundary appears. So we really observe that the two 
separate phase boundaries can be changed into one by the 
non-anomaly flavor-mixing induced by the mismatched vec- 
tor interactions. 

The above calculation for ôu = —20 MeV is further ex- 
tended to a fixed moderate coupling g; = 0.4g;. The phase 
diagrams for varied gý with g}; > g} are shown in Fig. 2, which 
is still analogous to Fig. 1. In contrast to Fig. 1, a stronger gy 
is required for the conversion of the two phase transitions into 
one due to the enlarged gj. Fig. 2 also shows that the chiral 
transition is first softened and then strengthened with g}. By 
comparison, the chiral transition is always weakened with the 
increase of g5. 

So for the weak isospin asymmetry, Figs. 1 and 2 show 
that the chiral transition separation can be removed by the 
mismatched vector interactions, even without the instanton 
induced flavor-mixing. Actually, all the two sets of phase di- 
agrams in Figs. | and 2 are quite similar to Figs. 2 in Ref. [23] 
obtained by changing the a. In this sense, the non-anomaly 
flavor-mixing due to the vector coupling difference plays the 
similar role as the KMT interaction. 

However, Figs. 1 and 2 indicate that gù must be much 
stronger than gý for turning the two chiral transitions into one: 
gy is at least twice as strong as g} to remove the separation. 
Of course, the fate of the separate chiral transitions depends 
on not only the vector coupling difference, but also the mag- 
nitudes of gy and gẹ. Here we do not show the results for 
gy > g5 with a fixed strong g; since in this case only crossover 
transition appears. 

On the contrary, we do not find the coincidence of the de- 
tached phase boundaries for gy < gẹ. In Fig. 3, we show 
the phase diagrams for du = —20 MeV with varied gi and 
fixed coupling g} = 0.2g,. We see that the two separate phase 
boundaries get farther rather than closer with the increase of 
|6g,| for g} < g5, which is quite different from what shown in 
Figs. 1 and 2. 

The reason can be traced back to Eqs. (14) and (15). First, 
according to Eq. (15), the |u;| is explicitly less than the |;| 
since the signs of uz and —2g}(p, — pa) in pi are different for 
g > 0. So for g; > gi with a fixed g5, increasing g}, implies 
not only the enhancement of the flavor-mixing but also the 
reduction of |y;|. This is why the two phase boundaries ap- 
proach each other with g}, as shown in Figs. 1 and 2. Second, 
near the left side of the right phase boundary, the p4 is remark- 
ably larger than the p, because of the significant suppression 
of the d quark mass; but around the left side of the left phase 
boundary, the difference between the pg and p, is relatively 
small. So for g5 > gi, the flavor-mixing term —(g* — g!)oq 
in uj, impacts the right phase boundary more significantly in 
contrast to what the corresponding term —(g} — g})pu in u} 
does on the left phase boundary, according to Eq. (14). This 
is why the right phase boundary moves more rapidly towards 
the higher u with gẹ in contrast to the left one, as shown in 
Fig. 3. 


$20511-5 


ZHANG Zhao and SU Hai-Peng 


60 
gy = 0.4g, 
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= ôu = —20MeV 
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0 1 fi 1 
350 400 450 
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Fig. 2. (Color online) The T-u phase diagrams for varied vector- 
isovector coupling g? at du = —20 MeV without the axial anomaly. 
The vector-isoscalar coupling is fixed as gf = 0.4 (relative to the 


scalar coupling gs). The solid line stands for the first-order chiral 
boundary. 


If g; or/and g} are strong enough, the first-order chiral tran- 
sition will change into crossover and it would be no critical 
point. Owning to vector interactions, it is possible that one 
of the two phase boundaries first disappears while the other 
one still remains with the change of vector interactions (In 


T (MeV) 


õu= -20MeV 


0 
300 350 400 450 500 
H (MeV) 


Fig. 3. (Color online) The T- phase diagrams for varied vector- 
isoscalar coupling g5 at du = —20 MeV. The vector-isovector cou- 
pling is fixed as g} = 0.2 (relative to the scalar coupling gs). The 
axial anomaly is ignored. All the lines stand for the first-order chiral 
boundaries. 


contrast, the two critical endpoints always appear at the same 
temperature in Refs. [10, 11]). Such a case is really observed 
in Fig. 3 for very strong vector interaction g} = 1.0g,. Actu- 
ally, our numerical study suggests that the emergence of only 
one critical endpoint via this manner does not require very 
strong vector interaction when the weak KMT interaction is 
included. 

Here only the weak isospin asymmetry is considered be- 
cause |uz| is small in heavy-ion collisions. On the other hand, 
the quark matter may appear in the core of neutron star. In 
such a case, the magnitude of the difference between the u 
and d quark chemical potentials may be as large as 100 MeV 
due to the constraint of charge neutrality*. For the strong 
isospin asymmetry with a relatively large |u;|, we find that 
the separation of the chiral transition can not be removed by 
the non-anomaly flavor-mixing without considering the axial 
anomaly. However, the similar phase diagram as Fig. 3 is still 
observed for a proper choice of gY and g%. 


V. DISCUSSION AND CONCLUSION 


We have studied the influence of vector interactions with 
different coupling constants in the isoscalar and isovector 
channels on the possible separation of the chiral transition un- 
der the isospin asymmetry in a two-flavor NJL model, where 
the U(1)4 symmetry is assumed to be restored effectively near 
the phase boundary. 

We first show that, besides the argument from the em- 
pirically different nucleon and vector-meson couplings [33], 


4 Note that in another case with |uz| > m,/2 and vanishing ju, the pion con- 
densation emerges at zero and low tmeperatures[39—4 1] 
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the one-gluon exchange type interaction can also give rise to 
unequal vector interactions with g5 > g, in the MFA when 
including the Fock contribution. On the other hand, by ex- 
tending the work [24] to finite uz, we obtain the quite different 
vector coupling difference with g5 < g} from the constraints 
of lattice chiral curvatures at zero/small quark chemical po- 
tentials. We demonstrate that, similar to the mass-mixing 
induced by the KMT interaction, the density-mixing of two 
flavors is produced owning to the mismatched vector interac- 
tions. 


The role of the non-anomaly flavor-mixing on the chi- 
ral phase transition is investigated under the condition with 
weak isospin asymmetry. We find that to convert the two 
separate chiral transitions into one, gy must be significantly 


[1] Philipsen O. The QCD equation of state from the lat- 

tice. Prog Part Nucl Phys, 2013, 70: 55-107. DOI: 

10.1016/j.ppnp.2012.09.003 

Stephanov A M. Non-Gaussian fluctuations near the QCD 

critical point. Phys Rev Lett, 2009, 102: 032301. DOI: 

10.1103/PhysRevLett.102.032301 

Athanasiou C, Rajagopal K and Stephanov M. Using higher 

moments of fluctuations and their ratios in the search for the 

QCD critical point. Phys Rev D, 2010 82: 074008. DOI: 

10.1103/PhysRevD.82.074008 

Stephanov A M. On the sign of kurtosis near the QCD critical 

point. Phys Rev Lett,2011, 107: 052301. DOI: 10.1103/Phys- 

RevLett.107.052301 

[5] Gupta S, Luo X, Mohanty B, et al. Scale for the phase diagram 

of quantum chromodynamics. Science, 2011, 332: 1525-1528. 
DOI: 10.1126/science. 1204621 

[6] Gavai V R and Gupta S. Lattice QCD predictions for shapes 

of event distributions along the freezeout curve in heavy- 

ion collisions. Phys Lett B, 2011, 696: 459-C463. DOI: 

10.1016/j.physletb.2011.01.006 

Mohanty B. STAR experiment results from the beam energy 

scan program at RHIC. J Phys G Nucl Partic, 2011, 38: 

124023. DOI: 10.1088/0954-3899/38/12/124023 

Mitchell T J. The RHIC beam energy scan program: Results 

from the PHENIX experiment. Nucl Phys A, 2013, 904-905: 

903c—906c. DOI: 10.1016/j.nuclphysa.2013.02.161 

Mohanty B. Exploring the QCD phase diagram through high 

energy nuclear collisions: An overview. The 8th Interna- 

tional Workshop on Critical Point and Onset of Deconfinement 

(CPOD 2013), Napa, California, USA, Mar.11—25, 2013. 

[10] Klein B, Toublan D and Verbaarschot M J J. The QCD phase 
diagram at nonzero temperature, baryon and isospin chemi- 
cal potentials in random matrix theory. Phys Rev D, 2003, 68: 
014009. DOI: 10.1103/PhysRevD.68.014009 

[11] Toublan D and Kogut B J. Isospin chemical potential and 
the QCD phase diagram at nonzero temperature and baryon 
chemical potential. Phys Lett B, 2003, 564: 212-216. DOI: 
10.1016/S0370-2693(03)00701-9 

[12] Kitazawa M, Koide T, Kunihiro T, et al. Chiral and color su- 
perconducting phase transitions with vector interaction in a 
simple model. Prog Theor Phys, 2002, 108: 929-951. DOI: 
10.1143/PTP. 108.929 

[13] Hatsuda T, Tachibana M, Yamamoto N, et al. New critical point 


[2 


faery 


[3 


= 


[4 


= 


[7 


(a; 


[8 


fate 


[9 


= 


Nucl. Sci. Tech. 2, S20511 (2015) 


stronger than gi without the axial anomaly. In this situation, 
the non-anomaly flavor-mixing induced by vector interactions 
impacts the separation of chiral phase transitions in the sim- 
ilar way as the anomaly one induced by instantons: the two 
detached phase boundaries get closer first and then coincide 
with the enhancement of the flavor-mixing. 

Note that recently the Polyakov-Loop extended NJL model 
has been extensively used to investigate the thermal and dense 
properties of QCD. We stress that introducing the Polyakov- 
Loop dynamics does not qualitatively change our main con- 
clusions. In addition, our study can be directly extended to 
the quark meson model of QCD by incorporating the quark- 
vector-meson couplings. Especially, it is interesting to inves- 
tigate the role of the non anomaly flavor-mixing on the possi- 
ble quark-menson transition in neutron star [42]. 


induced by the axial anomaly in dense QCD. Phys Rev Lett, 
2006, 97: 122001. DOI: 10.1103/PhysRevLett.97.122001 

[14] Zhang Z, Fukushima K and Kunihiro T. Number of the QCD 

critical points with neutral color superconductivity. Phys Rev 

D, 2009, 79: 014004. DOI: 10.1103/PhysRevD.79.014004 

[15] Zhang Z and Kunihiro T. Vector interaction, charge neutrality 

and multiple chiral critical point structures. Phys Rev D, 2009, 

80: 014015. DOI: 10.1103/PhysRevD.80.014015 

[16] Zhang Z and Kunihiro T. Roles of axial anomaly on neutral 

quark matter with color superconducting phase. Phys Rev D, 

2011, 83: 114003. DOI: 10.1103/PhysRevD.83.114003 

[17] Kunihiro T, Minami Y and Zhang Z. QCD critical points and 

their associated soft modes. Prog Theor Phys Suppl, 2010, 186: 

447-454. DOI: 10.1143/PTPS.186.447 

[18] Pisarski D R and Wilczek F. Remarks on the chiral phase transi- 

tion in chromodynamics. Phys Rev D,1984, 29: 338-341. DOI: 

10.1103/PhysRevD.29.338 

[19] °t Hooft G. Computation of the quantum effects due to a four- 

dimensional pseudoparticle. Phys Rev D, 1976, 14: 3432- 

3450. DOI: 10.1103/PhysRevD. 14.3432 

[20] °t Hooft G. How instantons solve the U(1) problem. Phys Rept, 

1986, 142: 357-387. DOI: 10.1016/0370-1573(86)90117-1 

[21] Kobayashi M and Maskawa T. Chiral symmetry and 7-X 

mixing. Prog Theor Phys, 1970, 44: 1422-1424. DOI: 

10.1143/PTP.44.1422 

[22] Bratovic M N, Hatsuda T and Weise W. Role of vector in- 
teraction and axial anomaly in the PNJL modeling of the 
QCD phase diagram. Phys Lett B, 2013, 719: 131-135. DOI: 
10.1016/j.physletb.2013.01.003 

[23] Frank M, Buballa M and Oertel M. Flavor mixing effects on 
the QCD phase diagram at nonvanishing isospin chemical po- 
tential: One or two phase transitions?. Phys Lett B, 2003, 562: 
221-226. DOI: 10.1016/S0370-2693(03)00607-5 

[24] Bazavov A, Bhattacharya T, Buchoff I M, et al. The chiral tran- 
sition and U(1), symmetry restoration from lattice QCD using 
domain wall fermions. Phys Rev D, 2012, 86: 094503. DOI: 
10.1103/PhysRevD.86.094503 

[25] Cossu G, Aoki S, Fukaya H, et al. Finite temperature study of 
the axial U(1) symmetry on the lattice with overlap fermion for- 
mulation. Phys Rev D, 2013, 87:114514. DOI: 10.1103/Phys- 
RevD.87.114514 

[26] Aoki S, Fukaya H and Taniguchi Y. Chiral symmetry restora- 
tion, eigenvalue density of Dirac operator and axial U(1) 


$20511-7 


ZHANG Zhao and SU Hai-Peng 


[27] 


[28] 


[29] 


[30] 


[31] 


[32] 


[33] 


[34] 


anomaly at finite temperature. Phys Rev D, 2012, 86: 114512. 
DOI: 10.1103/PhysRevD.86.114512 

Asakawa M and Yazaki K. Chiral restoration at finite den- 
sity and temperature. Nucl Phys A,1989, 504: 668-684. DOI: 
10.1016/0375-9474(89)90002-X 

Takizawa M, Kunihiro T and Kubodera K. Axial vector 
matrix elements Of the proton in the Nambu-Jona-Lasinio 
model. Phys Lett B, 1990, 237: 242-247. DOI: 10.1016/0370- 
2693(90)91436-F 

Klevansky P K. The Nambu-Jona-Lasinio model of quantum 
chromodynamics. Rev Mod Phys, 1992, 64: 649-708. DOI: 
10.1103/RevModPhys.64.649 

Zhang Z. Correction to the chiral magnetic effect from axial- 
vector interaction. Phys Rev D, 2012, 85: 114028. DOI: 
10.1103/PhysRevD.85.114028 

Cea P, Cosmai L, D’Elia M, et al. The critical line of two- 
flavor QCD at finite isospin or baryon densities from imagi- 
nary chemical potentials. Phys Rev D, 2012, 85: 094512. DOI: 
10.1103/PhysRevD.85.094512 

D’Elia M and Sanfilippo F. Thermodynamics of two flavor 
QCD from imaginary chemical potentials. Phys Rev D, 2009, 
80: 014502. DOI: 10.1103/PhysRevD.80.014502 

Sasaki C, Friman B and Redlich K. Quark number fluctuations 
in a chiral model at finite baryon chemical potential. Phys Rev 
D, 2007, 75: 054026. DOI: 10.1103/PhysRevD.75.054026 
Ferroni L and Koch V. Mean field approach to flavor sus- 
ceptibilities with a vector interaction. Phys Rev C, 2011, 83: 


[35] 


[36] 


[37] 


[38] 


[39] 


[40] 


[41] 


[42] 


$20511-8 


Nucl. Sci. Tech. 2, 20511 (2015) 


045205. DOI: 10.1103/PhysRevC.83.045205 

Tao C, Ma Y G, Zhang G Q, et al. Symmetry energy depen- 
dence of the pygmy and giant dipole resonances in an isospin 
dependent quantum molecular dynamics model. Nucl Sci Tech, 
2013, 24: 030502. 

Xie W J and Zhang F S. Probing the density dependence of the 
symmetry energy with central heavy ion collisions. Nucl Sci 
Tech, 2013, 24: 050502. 

Zhang Y X, Li Z X, Zhao K, et al. Constraints on symmetry 
energy and n/p effective mass splitting with improved quantum 
molecular dynamics model. Nucl Sci Tech, 2013, 24: 050503. 
Feng Z Q. Isospin splitting of nucleon effective mass and 
isospin effect in heavy-ion collisions. Nucl Sci Tech, 2013, 24: 
050504. 

Son T D and Stephanov A M. QCD at finite isospin den- 
sity. Phys Rev Lett, 2001, 86: 592-595. DOI: 10.1103/Phys- 
RevLett.86.592 

Zhang Z and Liu Y X. Coupling of pion condensate, chiral 
condensate and Polyakov loop in an extended NJL model. Phys 
Rev C, 2007, 75: 064910. DOI: 10.1103/PhysRevC.75.064910 
Zhang Z and Liu Y X. Two-flavor QCD phases and conden- 
sates at finite isospin chemical potential. Phys Rev C, 2007, 
75: 035201. DOI: 10.1103/PhysRevC.75.035201 

Shao G Y, Colonna M, Di Toro M, et al. Isoscalar-vector inter- 
action and its influence on the hadron-quark phase transition. 
Nucl Sci Tech, 2013, 24: 050523. 


